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Leading and subleading twist transverse momentum dependent parton distribution functions 
(TMDs) are studied in a quark model framework provided by the bag model. A complete set of 
^— ■>) . relations among different TMDs is derived, and the question is discussed how modeI-(in)dependent 

' such relations are. A connection of the pretzelosity distribution and quark orbital angular momen- 

I tum is derived. Numerical results are presented, and applications for phenomenology discussed. In 

C , particular, it is shown that in the valence-x region the bag model supports a Gaussian Ansatz for 

the transverse momentum dependence of TMDs. 
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p . I. INTRODUCTION 

^ TMDs are a generalization [H-Q of parton distribution functions (PDFs) promising to extend our knowledge of the 

] nucleon structure far beyond what we have learned from PDFs about the longitudinal momentum distributions of 
partons in the nucleon. In addition to the latter, TMDs carry also information on transverse parton momenta and 
spin-orbit correlations [5l-l35|. Here longitudinal and transverse refers to the hard momentum flow in the process, for 
^ , example, in deeply inelastic lepton nucleon scattering (DIS) the momentum of the virtual photon. 

' TMDs (and/or transverse momentum dependent fragmentation functions) enter the description of leading-twist 
\^ , obscrvablcs in deeply inelastic reactions [5|-[3l on which data are available like: semi- inclusive DIS (SIDIS) |36l - [50| . 
' DrcU-Yan process j5ll - [53j . or hadron production in e+e~ annihilations [5^ - l57| . 

] The interpretation of these data is not straight-forward though. In SIDIS one deals with convolutions of a priori 
7-H ■ unknown transverse momentum distributions in nucleon and fragmentation process, and in practice is forced to assume 
■ models for transverse parton momenta such as the Gaussian Ansatz [58l - l67| . In the case of subleading twist observables, 
. one moreover faces the problem that several twist-3 TMDs and fragmentation functions enter the description of one 
observable [68l - l77j (we recall that presently factorization is not proven for subleading- twist observables [TSj). 

In this situation information from models [t^ - Iq^ ] is valuable for several reasons. Models can be used for direct 
estimates of obscrvablcs, though it is difficult to reliably apply the results, typically obtained at low hadronic scales, 
to experimentally relevant energies (86j . Another aspect concerns relations among TMDs observed in models (tS * 



- ' — \ 

X 



■ Such relations, especially when supported by several models, could be helpful — at least for qualitative interpretations 
of first data. Furthermore, model results allow to test assumptions made in literature, such as the Gaussian Ansatz 
for transverse momentum distributions or certain approximations [osl - llOl | . 



In addition to such practical applications model studies are of interest also because they provide important insights 
into non-perturbative properties of TMDs. In this context the probably most interesting recent observation in models 
concerns the pretzelosity distribution function, which in some quark models is related to the difference of the hclicity 
and transversity distributions (s^l and, so far, in one model to quark orbital momentum [sj] which is, to best of our 
knowledge, the first 'rigorous' connection of a TMD and quark orbital angular momentum in a model. 

The purpose of this work is to study TMDs in the framework of the MIT bag model. We compute in this model 
all leading- and subleading- twist, time-reversal (T-) even TMDs in Sec.|lTl and address then in Sec. IIIII questions like: 
how do relations among TMDs arise in a quark model? How many such relations are there in a model? To which 
extent may one expect such relations to be realized in nature? In Sec. IIVI we establish a connection of pretzelosity 
and quark orbital angular momentum in the bag model. In SeclVlwe present and discuss the numerical results, using 
them, among others, for 'testing' the Gaussian Ansatz or Wandzura-Wilczck-type approximations |98l4lOll . Finally, 
in Sec. IVII we present our conclusions. Some of the results presented here were shown in the proceeding [871 . 

For convenience and in order to make this work self-contained, in the remainder of this Introduction we include 
general definitions of TMDs, and introduce relevant notation. 
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A. General definitions of TMDs 



Hard processes sensitive to parton transverse momenta like SIDIS are described in terms of light-front correlators 



(27r)3 



g.p. (iV(P,5)|7^,(0) W(0, z; path)^,(z)|7V(P,5)) 



(1) 
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We use lieht-cone coor dinatcs a± = (a° ±a^)/y/2. In SIDIS the singlcd-out 3-direction is along the momentum 
of the hard virtual photon, and transverse vectors like pV are perpendicular to it. The path of the symbolically 
indicated Wilson-link depends on the process j23l - l25| . In the nucleon rest frame the polarization vector is given by 

S = (0, St, Sl) with 5^1,-^51, = 1. 

The information content of the correlator ([T|) is summarized by eight leading- twist TMDs [3], that can be projected 
out from the correlator (P) as follows (color online: red: T-odd, blue: T-even) 
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where the space-indices j, k refer to the plane transverse with respect to the light-cone and e^^ = —e^^ ~ 1 and zero 
else. Integrating out transverse momenta in the correlator ^ leads to the ' usual ' par ton distributions known from 
collinear kinematics i°'{x) ~ J (PpT j'^ix^pp) with j = fi, <7i, hi, e, ,gT, h^ jl02l Il03| . Dirac-structures other than 
that in Eqs. (pi-[TU| lead to subsubleading-twist terms p^. 

For convenience we introduce for a generic TMD j'^(x,k±) the '(unintegrated) transverse (l)-moments' defined as 



j(i)«(x, k^) = fix, fcx) , /i)«(x) 



2APn 
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Moreover, we shall also make use of the '(l/2)-moments' defined for a generic TMD as 



(11) 



(12) 
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II. TMDS IN THE BAG MODEL 



In the MIT bag model, the quark field has the following general form 



*a(^,i)= X! N{nK){ba{nKm)il;nKj^{x,t) + dl^{nKm)i;_ri-K.jTnix,t)} , (13) 

n>0,K;=±l,m=±l/2 

where &J, and create quark and anti-quark exeitations in the bag with the wave functions 

^„.-i.i™(-^,i) = ^f . ) (14) 

For the lowest mode, we have n = 1, k = —1, and ~ 2.04 denoted as w = tui^-i in the following. In the above 

equation, a is the 2x2 Pauli matrix, Xm the Pauli spinor, Rq the bag radius, x = x/\x\, and ji, are spherical Bessel 
functions. Taking the Fourier transformation, we have the momentum space wave function for the lowest mode, 

^^{k) = ^V4^NRl( ^ S.!°JS?0 ' (15) 



a ■ k ti(k)x 

where k = k/k with k ~ \k\ and the normalization factor N is 



s \ 1/2 



T7 ^ ^ • (16) 

The two functions ti, i = 0,1 are defined as 

1 

U{k) = / u'^ duji{ukRo)ji{uu) . (17) 



From the above equations, we see that the bag model wave function Eq. (|15p contains both S and P wave components. 
Especially, to represents the S'-wave component, whereas ti represents the P-wave component of the proton wave 
functions. 

With the above wave functions, we can calculate all quark TMDs. For convenience we define the constant A, which 
will be common to all TMDs, and the momenta k^ and k as 

^ = -2? N,f2 ' k = Jkl + kl , k, xMn - Lo/Ro , fc. = , M)v = ^ , (18) 



where Mn is the proton mass, and the bag radius is fixed such that RqMn = 4a;. Moreover, we assume SU{6) 
spin-flavor symmetry of the proton wave function, such that spin-independent TMDs of definite flavor are given in 
terms of respective 'flavor- less' expressions multiplied by a 'flavor factor' Nq, and spin-dependent TMDs of definite 
flavor follow from multiplying the respective 'flavor-less' expressions by a 'spin-flavor factor' Pq with 

iV„ = 2 , Nd = l, Pu^^, = "J ■ 

We recall that in the quark model formulated for a general (odd) number of colors Nc, these flavor factors are given 
by Nu = (iVc + l)/2 and = {N, - l)/2 while P,, = {N, + 5)/6 and P^ = {-N, + l)/6 

We mention that the MIT bag model gives rise also to antiquark distributions, but to unphysical ones, since 
fi{x) < 0, which violates positivity. The TMDs receive no n- vanishing support also from the regions \x\ > 1. Though 
non-physical these contributions must be included when evaluating sum rules like J dx f1{x) = Ng or the momentum 
sum rule J dxxff{x) = 1, i.e. sum rules are satisfled only when integrating over the whole x-axis. 

In literature it was discussed how to deal with these caveats, see for example |108| . In this work, we limit ourselves 
to the discussion of quark TMDs at < a; < 1, which should not be confused with 'valence distributions', for example 
/ivai(^) ~ /i (^) ~ /i (^)- When discussing sum rules, however, integration over the whole cc-axis is implied. 

Since there are no explicit gluon degrees of freedom, T-odd TMDs vanish in this model (9l| . In principle, one can 
simulate the effect of the gauge link, which is crucial in QCD for T-odd effects [22| - [2^ . for example by introducing 
'one-gluon-exchange' [9l|, |92| or invoking instanton effects [Q^I • In this work we shall not consider such extensions of 
the bag model, and restrict ourselves to the description of T-even distributions. 



A. Results for TMDs in the bag model 

In the notation introduced above, the results for the T-even leading twist TMDs are given by 
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and for the subleading twist TMDs we obtain 
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In the following Seetions we shall discuss these results in detail. 
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III. EQUALITIES AND INEQUALITIES AMONG TMDS 

In QCD all TMDs arc independent functions. But in quark models, due to absence of gauge field degrees of freedom, 
certain relations among different TMDs appear which must be satisfied in any consistent relativistic quark model. 
We discuss these "model-independent" quark-model relations in Sec. IIII Al Of course, depending on a quark model 
further relations may appear, and the bag model results (|20j|33p provide a nice illustration why this happens which is 
demonstrated in Sees. IIII Bl and HTl Cl in detail. In Sec. IIII Dj we compare to results from other models. This comparison 
helps to establish to which extent which relations might be expected to be useful in nature. Sec. IIII El is devoted to 
the discussion of one particular relation. 



A. Relations valid in all quark models 

Ce rtain relations among TMDs must be valid in any quark model of the nucleon lacking giuon degrees of freedom 
[lOCf . In such "no-gluon models" the absence of the Wilson- link implies that in the general Lorentz-decomposition of 
the unintegrated quark-correlator certain amplitudes do not appear, namely the i3i-amplitudes (i = 1, 2, . . . 20) in 
the notation of [2^1 . This gives rise to the following relations [13, llOCll | 

grix) W g^^:c) + ^gi^'\x) , (34) 
hL{x) h,{x) - -^hfl^'\x) , (35) 



da; 

hrix) ™ -^hfi'\x), (36) 



9iix) + ^9^^'\x) 0, (37) 

hT{x,pT) - h^{x,pT) ^= h^L{x,pT), (38) 

which must hold in any consistent relativistic quark model. These so-called "Lorentz-invariance relations" (LIRs) 
are not valid in models with gauge field degrees of freedom [o^ and in QCD [9^. The applications of LIRs in 
phenomenology were discussed in |99l . Il00l |. There it was also shown, by exploring QCD equations of motion, that 
some LIRs hold in an approximation consisting of the neglect of quark-gluon-quark-correlator and current quark mass 
terms. Whether suc h an approxim ation is is justified in nature is, of course, a different question. For discussions of 
specific cases see [98l4lOll Il09l4ll3 |. For quark model calculations, the practical application of the relations 



is immediate: they provide a valuable cross check for the numerical results. 

In App. |X] we provide analytical proofs that the LIRs ([MH55| are satisfied in the bag model. We also checked that 
the numerical results satisfy the LIRs, which provides a welcome cross-check for the numerical calculation. 
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B. Linear relations in bag model 



In the bag model, there are 9 linear relations among the 14 (twist-2 and 3) T-evcn TMDs, which can be written as 
follows 



e«(x, fcx) + hl{x, k^) = 2.g|,(x, 
P'?/^9(a;,fci) = 4^(.T,fcx) 
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where the 'dilution factor' is defined as 
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In the relations (|45l |46)) some TMDs need to be multiplied by the model-independent factor fc^/(2Af^), which is 
a 'legitimate linear operation' in our context (the meaning of that will be explained shortly). The '(unintegrated) 
transverse moments' are defined in Eq. ()lip . 

Why are there 9 linear relations? In fact, naively, one could have expected even more relations, since all TMDs are 
expressed in terms of only two functions, Iq and ti representing the contributions from the S and P-wave components 
of the proton wave function, Eqs. (|15l I17p . However, having linear relations in mind, the combinations t^, io^ii 
are to be considered as independent structures. But there are more independent structures than that. By inspecting 
Eqs. (|20j|33p we see that the actual number of linearly independent structures in the TMDs is 5, namely 
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(49) 



where we show respectively to which TMDs the different structures contribute. We observe that in some sense there 
are two 'disconnected subspaces': one is due to the structures I, II, III, and the other due to the structures IV, V. 

The structures II and IV, fc^io^i and t^ti, are linearly independent, as there is no way of relating one with the other 
in a model-independent way. Indeed, in order to do this, one should multiply a TMD by a factor including kz which 
explicitly depends on parameters of the bag model, as is evident from Eq. (|18p . and we discard such a manipulation 
as a model-dependent operation. For the same reason the structures in III and V are linearly independent. 

But k'^t'l and in point III are linearly dependent: if we multiply (actually in relevant expressions M]qt\ appears) 
by the model-independent factor k^/Mf^ and add fc^ t\ we obtain just t\ which happens in Eqs. (|^ Clearly, the 
multiplication of TMDs by fc5^/(2M^) is a model-independent manipulation leading to transverse moments in (jlip . 

To summarize, there are 5 linearly independent structures in the bag model and 14 T-even TMDs. This implies 
9 linear relations, and Eqs. (|39H47|) represent one way of writing these relations. These findings mean that one can 
choose, in the bag model, a basis of 5 linearly independent TMDs, and construct the other TMDs from this basis. 
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C. Non-linear relations in the bag model 



In Eq. (|49| we found that the TMDs in the bag model form two independent 'subspaces.' Let us summarize: 

suhspace A {I, II, III): ff gl hi hf^ e« 5* /i« 5^'' => relations (I2S1 SOI 113 Hi Hi) , (50) 



twist 2 twist 3 

subspace B (IV, V): g^-^ h^l f^" 

twist 2 twist 3 



relations gH HI HI ST]) . 



In other words, there is no linear relation which would transform TMDs from subspaee A into TMDs in subspace B. 
However, there are non-linear relations which can do that, for example. 
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The results are presented such that on the left-(right-)hand-sides (L(R)IIS) only TMDs from subspace A (B) appear. 
The Eqs. ((5TI [52|) are independent in the sense that it is impossible to convert one into the other upon use of the 



linear relations ([5^^T7|) . In order to convince oneself of that, notice that on RHS of ([51]) h^^ £ subspace B appears 
from which alone one cannot construct the TMDs £ subspace B on RHS of (|52l) : from h^^ one obtains g^^ via 
but not gj^''- However, for example, g^^ and gj^"^ span a basis which allows to construct all TMDs € subspace B. 
Similarly one finds that from the over-complete set of TMDs on LHS of ([FTl [5^ all TMDs e subspace A follow. To 
summarize, the non-linear relations (j51[l52[) are independent, and these are the only independent non-linear relations. 

Of course, upon the use of the linear relations (|5Wff7)) one could generate further non-linear relations. One advantage 
of the presentation ([5T1 [5^ is that they connect only chirally even, or only chirally odd TMDs. 

With the 9 linear relations ([551 - H71) and the 2 non-linear relations (|51[ [5^ we find altogether 11 relations among 
14 TMDs in the bag model. This may refiect that eventually all TMDs can be traced back to the free structures 
proportional to Iq, toti, t^. 

It should be noticed that in Sees. IIII Bl and IIII CI we permitted only manipulations of the kind: adding TMDs, 
multiplying them, and forming (l)-moments. If one includes differentiation of TMDs one obtains additional relations, 
see Sec. IIII Al and App. 
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D. How general are quark model relations among TMDs? 



The deeper reason, why in the bag model relations among TMDs appear, is ultimately related to Melosh r otat ions 



which connect longitudinal and transverse nucleon and quark polarization states in a Lorentz-invariant way |114 |. 

An important issue, when observing relations among TMDs in a model, concerns their presumed validity beyond 
that particular model framework. For that it is instructive to compare first to other models. In fact, some of the 
relations (|39H47p were discussed previously in literature in various models. Let us review briefly. 



• Eq. ((39|) : i ts k i -integrated version was discussed in bag model in jl03j and jll5l . lll6t and in light-cone constituent 
models in [llTj . The unintegrated version was discussed in bag and light-cone constituent models {sOi. iSli]. 

• Eq. (I40p : its integrated version was observed in the bag model previously in [ll5j | . 

• Eq. (|42p : was first observed in the spectator model of [t^ and later also in light-cone constituent models [8l| 
and the covariant parton model of Ref. [Hj . 

• Eq. (H^ : was found in the spectator model of Ref. [t^ . 

• Eq. (gS]) : was first observed in the bag ^ . It is valid also in the spectator [t^ , light-cone constituent [8l| , and 
covariant parton f83| models. 

• Eqs. (|411 1151 H51 ITT]) : are new in the sense of not having been mentioned previously in literature. But the latter 
3 are satisfied by the spectator model results from [79| . 

• The non-linear relation ([?T|) . which connects all T-even, chirally-odd leading- twist TMDs was observed in the 
covariant parton model approach (s^ . Eq. (j52p was not discussed so far in literature. 

The detailed comparison, in which models these relations hold and in which they are violated, gives some insight into 
the question to which extent these relations are model-dependent. 

Let us discuss first Eqs. ([5M1T|) . which include the 'dilution factor' (|15)) and connect polarized and unpolarized 
TMDs. For these relations SU(6)-spin-flavor symmetry is necessary, but not sufficient. In fact, this type of relations 
holds only in 'simplest models' such as the bag model version used here or light-cone constituent models [10, [8l| . 
What these models have in common is that the nucleon wave-function is constructed from 'flavor-blind' quark wave- 
functions multiplied by appropriate spin-flavor factors in Eq. (|19p. The SU(6) symmetry, however, does not need to 
be realized in a model that simply. For example, the spectator model of [79| is SU(6) symmetric. But it does not 
support p9l - |4T|) which are spoiled by the different masses of the (scalar and axial- vector) spectator diquark systems. 
Interestingly, it is possible to recover these relations in [t^ in the limiting case of the scalar and axial-vector diquark 
masses becoming equal (justified in large- A^c limit). We mention that ([551 HO]) also are not supported in the covariant 
parton model approach of (s^ . However, also in that approach it is possible to 'restore' these relations by introducing 
additional, restrictive assumptions, see [83 for a detailed discussion. We conclude that the relations (|39fl4T|) require 
strong model assumptions. It is difficult to estimate to which extent such relations could be useful approximations in 
nature, though they could hold in the valence-x region with an accuracy of (20-30) %, see the comparison of similar 
SU(6) predictions and data for the polarized neutron and proton structure functions in [86j . 

From the point of view of model dependence, it is 'safer' (sot to compare relations which include only polarized or 
only unpolarized TMDs. We know no example for the latter, however, the relations (|42||T7|) are of the former type. 
It is gratifying to observe that these relations are satisfied not only in the bag model, this work and [s^, but also 
in the spectator model version of Ref. [t^. The relations among the leading twist TMDs, Eqs. ([^ H5|) . hold also in 
light-cone constituent [8l|, and covariant parton [s^l models. They are also valid in the non-relativistic model [s^ . 
Though they do not prove anything, these observations indicate that such relations could be valid in a wider class of 
quark models. 

Of course, quark model relations among TMDs have limitations, even in quark models. In [s^ various versions 
of spectator models were used, and in some versions the relations were not supported (|42i I45p . It is instructive to 
observe that also the quark-target model [s^ not supported the relations (|^ In fact, the inclusion of gauge 
fields brings us a step closer to QCD as compared to quark models, at least in what concerns the involved degrees of 
freedom. Finally, in QCD none of such relations is valid, and all TMDs arc independent structures. 

It would be interesting to 'test' such quark model relations in other models, lattice QCD, and in experiment. The 
latter, in fact, provide a test for the usefulness of quark models themselves — or, more precisely, their application to 
TMDs. 
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E. A special linear relation among collinear functions 

By taking linear combinations of P9H47[) one can obtain many more linear relations. It is worth to discuss in some 
more detail one particularly interesting relation, which can be obtained in this way. By eliminating the transverse 
moment of the pretzelosity distribution from Eqs. (|^ HS]) . and integrating over transverse momenta, we obtain 

glix) - hUx) = gU^) - hlix) . (53) 

This relation holds also in its unintegrated form. There are several reasons, why this relation is interesting. 

First, it involves only collinear parton distribution functions, which is the only relation of such type in bag model. 
(Actually (p9l |40| are also of such type, but they include the 'dilution factor' and are supported only in models with 
simplest spin-flavor structures, see Sec. IIII Dl ) The QCD evolution equation for all these functions are known, and 
they are different, which shows the limitation of this relation: even if for some reason (j53p was valid in QCD at a 
certain renormalization scale fxo, it would break down at any other scale fi ^ fiQ. This is by no means surprising, and 
we expect such limitations for all model relations. 

Second, for the first Mellin moment this relation is valid model-independently. Hereby we strictly speaking presume 
the validity of the Burkardt-Cottingham sum rule, which is equivalent to the statement J dx g!^{x) = J dx gKx), and 
an analog sum rule for h\{x) and h\{x). In QCD there are doubts espec ially concerning the validity of the Burkardt- 
Cottingham sum rule. However, it is valid in many models such as bag |l03l | or chiral quark soliton model jll8j |. 

Third, though it certainly is not exact in QCD, it would be interesting to learn whether (j53p is satisfied in nature 
approximately. Also this relation can be tested on the lattice, especially for low Mellin moments and in the flavour 
non-singlet case. Lattice QCD calculations for Mellin moments of g'^{x) were reported in |112| |. 

Forth, the relation (|53p can be tested in models where collinear parton distribution functions were studi ed. Some 
results can be found in literature. For example, c alcul ations of parton distribution functions in bag models jl03l Ill5j | 
support this relation (the bag model version of jl03l | coincides with the one used here) . Moreover, the spectator 
model [Z^ supports this relation: it is equivalent to g\{x) = ^h\{x) in the notation of [t^], and also the unintegrated 
version of (1531) is valid there. One counter-example is known though: the chiral quark-soliton model does not support 
this relation |ll8l . Il26| . This observation could provide a hint in which models ((53)) is valid. The models where 
(|53p holds include only the components in the nucleon wave-function with the quark orbital angular momenta up to 
i = 0, 1, 2 at most. The chiral quark soliton model, which does not support ([53]). contains all quark angular momenta 
L = 0, 1, 2, 3, 4, . . . but this point deserves further investigation. 

Fifth, an important aspect of model relations is that they inspire interpretations. The relation (|53p means that the 
difference between and /i^ is t o the same extent a 'measure of relativistic effects in the nucleon' as the difference 
between helicity and transversity |103| . Both these differences are related to the transverse moment of pretzelosity, 
see Eqs. (|i5lH5| and [sot. 

F. Inequalities 



Finally, we discuss inequalities among leading twist TMDs, which are valid in QCD and all models [21 1 



/f(x,A:i)>0, \gl{x,k^)\<fl{x,k^), \h\{x,k^_)\ < fl{x,k^) , (54) 

\hl{x,k^)\ < i(^/f(x,fcx)+5?(a^,fc±)) , (55) 

\ht^{x.k^)\ < i(^/f(x,fcx)-5?(2^,fc±)) , (56) 

gt^'^^x^k^f + ft^'^^x^k^Y < ^ (jl{x,k^f ~ glixMf^ , (57) 

h^i'^%x, k^f + ht^'^\x, k^f < (^nix, k^f - gfix, k^)'^ , (58) 

where we have to keep in mind that in the present quark model framework the inequalities simplify, due the absence 
of the T-odd TMDs /{^ and h^'^ . In App. |B] we demonstrate explicitly that the bag model expressions for the 
quark TMDs satisfy ([SlHSS]) . It is interesting to remark that for the nucleon, except for f1{x,k±) > 0, all the other 
inequalities in ([51H55)) are 'true' (i.e. never saturated) inequalities, see App. iBjfor a detailed discussion. 
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IV. PRETZELOSITY AND QUARK ORBITAL ANGULAR MOMENTUM 



In quark models, in contrast to gauge theories, one may unambiguously define the quark orbital angular momentum 
operator as = ipqe'^^^f^p^ipq where for clarity the 'hat' indicates a quantum operator. This definition follows (in the 
absence of gauge fields) uniquely, for instance, from identifying that part of the generator of rotations not associated 
with the intrinsic quark spin. For the following it will be convenient to introduce a 'non-local version' of this operator, 
by defining (keep in mind that we work in non-gauge theory) 



L;(0,z) = V^,(0)e^'='fVV',(^). 

In the bag model it is convenient to work in momentum space where 
analogy to Eq. ([T]) the following quantity 



(59) 

and — ■ Next let us define, in 



L^ix^pr) 



dz- 



(27r)3 



{N{P,S')\LliO,z)\N{P,S^)) 



(60) 



z+=0,p+=xP+ 



In order to find a connection to TMDs we must consider a longitudinally polarized nucleon, choosing the polarization 
vector as iS* = (0, 0, 1) for definiteness, and we must focus on the j = 3 component in (|60p . i.e. on the component of 
the angular momentum operator along the light-cone space-direction. This is because the transverse momenta pt of 
the quarks generate orbital angular momentum which is oriented perpendicular to px (and to the transverse position 
of quarks inside the nucleon, which can be quantified rigorously in the impact parameter space in terms of generalized 
parton distribution functions, but we do not need this notion here). 

In a quark model, where the ambiguities of gauge field theories are absent, the partonic interpretation of (|60)) is 
the following. For example, in a longitudinally polarized nucleon Lq{x,pT)d^PTdx tells how much the orbital angular 
momentum of a quark of fiavour q, which carries the longitudinal momentum fraction x and the transverse momentum 
Pt = \pt\i contributes to the nucleon spin. (In QCD such an interpretation for the light-cone plus-component L'^ 
would also be possible, in an appropriately fixed gauge.) 

Evaluating the expression (pO]) in the bag model we obtain 



Llix,pT) = {-l)ht^'^''{x,pT) 



(61) 



In order to demonstrate the consistency of this result we compute the contribution to the total angular momentum 
of the nucleon due to flavour q. is composed of contributions from intrinsic quark spin, ~ ^ J dxgKx), and 
quark orbital angular momentum Lq= J dx J d^pTLq(x,pT)- We obtain 



2J„ = 



dx I d^k 
A 



« M. 



N 



P. 



gl{x,k^) - 2h^^^'^\x,ki_) 



d:'k 



/j2 

if, + 2kz toh + {2kl - 1 + 2 -i) 



/fc2 



(62) 



where we first substituted x — >■ fc^ = xMj^ — uj/R (recalling that a;-integration is carried over entire a;-axis. Sec. Ill|. 
then used that under the integral over d^fc for symmetry reasons k^ = k^/k drops out while fc^ — ?> i fc^ and fc^ § 
and finally observed the same integral which appears in the normalization of the unpolarized distribution. 

Eq. (|62p is the correct SU(6) quark model result for the contributions of various fiavours to the total nucleon spin. 
Notice that j3 J| = i. This confir ms that the connection of pretzelosity and the quark orbital angular momentum 
content of the nucleon is consistent. Thus, our results, supported by the light-cone SU(6) quark-diquark model [s^ . 
suggest that 

Ll = (-1) J dxhf^^^^ix) . (63) 

It is important to observe that the relation of pretzelosity and orbital angular momentum, Eqs. (|6ip and (|63p . is at the 
level of matrix elements of operators, and there is no a priori operator identity which would make such a connection. 
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V. NUMERICAL RESULTS 



In this Section wc discuss the numerical rcsuhs for the TMDs. In Sec. IV Al we present the results for the integrated 
TMDs as functions of x. In Sec. IVBl we focus on their /c^-behaviour. Finally, in Sec. IV Cl we investigate the question 
whether the bag model results support the so-called Wandzura-Wilczek(-type) approximations. 



A. Results for the integrated TMDs 



As the flavour dependence governed by the spin-flavour SU(6) symmetry is trivial, we will show only results for the 
M-flavor. For unpolarized TMDs the d-quark distributions are factor 2 smaller than the u-quark distributions. In the 
case of the polarized TMDs, the d-quark distributions are factor 4 smaller and have opposite sign compared to the 
M-quark distributions, according to Eq. (|T9| . All results discussed below refer to the low scale of the bag model. 

Let us start the discussion of the numerical results with unpolarized TMDs. Fig. [1^ shows the twist-2 unpolarized 
distribution function fi{x), and the subleading twist functions /-^"{x), e"{x). Only /^'(x) = / d'^k±f-^'^{x,k±) is 
new in this figure. The remarkable observation is that /-""'(x) is rather large, even larger than fi{x) for x < 0.7. 
However, one has to keep in mind that there are no positivity bounds for twist-3 TMDs. Moreover, it is ^f-^'^{x, k±) 
which enters in cross sections, and typi cally «C Mpf, which eventually guarantees positivity of cross sections . We 
remark that QCD equations of moti on |ll4l | imply a 5(x)-contribution_in effcc), which is found in some jll9l Il20| but 
not all effective approaches [tqI . I121| . including the bag model, see (l03l Ill5l | and Fig. [T^. 

Fig. [TJd shows the polarized functions gj^'^{x), g^^{x), h\{x), gi{x). The TMDs h^^{x) and hj^£{x) are simply 
related to the shown TMDs according to h^^{x) = —g^°'{x) and hj^£{x) = —g^rji^x), such that the results for these 
TMDs do not need to be shown. Also the results for d-quark distributions are not shown, as explained above. The 
results for the TMDs 5^"(a;), g^rp{x), h^£{x), are new, and it is interesting to observe that they are rather 

sizable, but again there are no positivity constraints on these objects. 

Fig. [1]; shows the polarized functions gj^^{x), (7^"(x), h'^{x), g^{x). The TMD h'^{x) is not shown, being related 
to gj;'^{x) as h'^{x) = —gj^'^{x). We see that hj^'^ is rather sizable, it is even bigger than /^(x) (the same scale is used 
in Figs. [TJi-c). Again there is no positivity constraint for this TMD, which would object this. 

The large size of the integrated twist-2 TMDs g^rji[x), h-^l{x), h-^'f:{x) can be understood qualitatively in the non- 
relativistic limit which was formulated for an arbitrary number of colours Nc in js^ . and can eventually be traced 
back the convention of using the nucleon mass in order to compensate the dimension of the k±^ factor (s) in the 
decomposition of the correlators in Eq. ([2l-[T0|). 

It is interesting to notice that g^''{x) = —h'ip{x) are the only TMDs in the bag model which have a zero in the 
valence-x region. This observation is actually not surprising but a consequence of the fact that the LIRs (PSIIST]) hold, 

and = —h'^^^'^ix) have extrema in the valence-x region. 



f\fl,e 



(a) 




8t' Sn' ^1' Si 



(b) 



, J- J- 



gj (x)=-hjj(x) 

J-M, , i-Lu, . 



bag model 




(c) 




0.2 0.4 0.6 0.8 



0.2 0.4 0.6 0.8 



0.2 0.4 0.6 0.8 



FIG. 1: (a) The unpolarized functions /^"(s), fi{x), e"(a;) vs. x from the bag model at the low scale. The d-quark distributions 
are factor two smaller compared to the unpolarized w-quark distributions according to the SU(6)-flavour factors in Eqs. P9II27|I . 

(b) The polarized functions g^'^{x) = —h^^{x), Qirix) = —h^£{x), hi{x), gi{x) vs. x. The d-quark distributions are factor 
four smaller and have opposite sign compared to the u-quark distributions according to the SU(6)-flavour factors in Eqs. (|19II27|I . 

(c) The polarized functions /iy"(a;), gj^^ix) = —h1j,(x), gj-ix), hf^ix) vs. x. The d-quark functions are as in (b). 
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B. Transverse momentum dependence 



In the context of TMDs the most interesting aspect is, of course, their transverse momentum dependence. In 
principle, all information is contained in the two-dimensional functions j{x, kj_) for a generic TMD, but here we shall 
content ourselves to discuss 'one- or zero-dimensional' projections of that information. 

The first point we address is: what are the typical transverse momenta of unpolarized quarks in the bag TMDs? 
For that we define for a generic TMD j'^{x, k±) the following quantities 

^^^^ ^ JdxJ^k^k^lixM} ^ ,^,,_ JdxJd'k^ klfjx^k^) ^ ^g^^ 
J da; / d^fcj^ fc^) ' J dx J d^fcj_ j(a;, fc^) 

Due to the simple spin flavor structure of the MIT bag model the {pr) and {p^) are flavor-independent for all TMDs. 

The flrst observation is that depending on the TMD {px) and {p^) in Eq. ([M)) may not exist in the bag model, 
because the momentum-space wave-function components ti(k), Eq. (|17p. do not vanish sufficiently fast at large k. 
This is the case especially for f1{x, k±). 

For the same reason also the (l)-moment /}^^*(a;) does not exist. However, the (l/2)-moment f[^^^^'^{x) defined 
according to (jl2p exists, and can be used to introduce an x-dependent average transverse momentum {pt{x)) as 

.(1/2)9. X 

{pt{x)) ^ 2Mn ^\ ^ . (65) 
J I \^) 

Fig. [2^ shows the result for f^^'^'^'' [x). (The divergence of {jpt) from (|64|) emerges when one tries to integrate 
over x, recalling that this integration extends to the entire x-axis, see Sec. Ull) 

Now the (l)-moment f^^'^{x) is divergent, but its derivative with respect to x exists, see the dotted line in Fig.[5]3. 
Hereby it is understood that the (l)-moment is computed with a finite cutoff Acut ^ Mm, then the derivative is 
taken, and only then the limit A^ut — >■ oo is performed. 

By integrating the well-defined ■^fi^^'^{x) we can compute a regularized (l)-moment {x)reg- The result depends 
on some arbitrary integration constant, which we fix such that the (l)-moment vanishes at a; = 1. This choice is 
reasonable but not unique, if we recall that in the MIT bag model TMDs in general have a non-zero (though small) 
support for |x| > 1, see Sec. [Hi Our main conclusions in this respect, to be presented below in this Section, depend 
weakly on the chosen value of the integration constant, provided reasonable choices are made (such as, for example, 

f[^^%x)rea = f[^^^^%x) at X = 1). The result for /f ^''(l)reg defined in this way is shown as solid line in Fig. [2b- 

With fi^'^'^ ix)reg wc arc in the position to define an x-dependent average transverse momentum square {p'^{x)) as 

(p|(x)) = 2M^ . (66) 




0.5 



-0.5 



(b) 




<Pj(x)> 



(c) 



0.5 
0.4 
0.3 
0.2 
0.1 




[-4<Pj{^)>\ 




bag model 



0.2 0.4 0.6 0.8 X 



0.2 0.4 0.6 0.8 X 







0.2 0.4 0.6 



FIG. 2: For the unpolarized TMD fi{x,k±) (a) the (l/2)-moment defined in Eq. (f65|) . (b) the derivative of the (l)-moment 
and the regularized (l)-moment as discussed in the text, and (c) {pt{x)) in comparison to {Tv{p^{x))/Ay^^. In the Gauss-model 
the two quantities would be equal. (The dotted marks the value {pt{x)) — 0.25 GeV.) 
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Fig-llt shows {pt{x)) as solid line. We observe that in the valence- a: region at the low hadronic scale {pt{x)) very 
weakly depends on x. Numerically we find 

{pt{x)) w 0.25 GcV for 0.2 < x < 0.5. (67) 

(The (pt(x)) — 0.25 GcV is marked as dotted line in Fig.[2t.) This is similar to results from the light-cone constituent 
model |8l| which also refer to a very low hadronic scale. In fact, keeping in mind the p^-broadening effects due to 
gluon radiation with increasing normalization scale [ll| . this is a reasonable result at a low scale. (We remark that 
in parton model approaches one finds comparably low values for {pt{x)) (albeit there the results refer to high scales) 
[12a [m models.) 

In phenomenology at high scales, however, larger values arc required [58l - [60| . For example, the interpretations of 
SIDIS data from EMC [IJ or HERMES [H require 

, , ^, ^r0.64GeV from EMC data in m, , . 

{Pt [x))gi,uss ~ j 0.56 GcV from HERMES data in [H, ^ ' 

where the index "Gauss" indicates that the Gaussian model has been assumed in these studies. The Gaussian model 
means that f^{x,pT) = fl{x) cxp(-p|/(p|(a;))Gauss)/(7r(py (x))Gauss)- The width (p|(x))Gauss could be a function 
of x, but in practice it is often assumed to be a constant. Such an Ansatz works with sufficient precision for many 
practical applications in phenomenological studies [58l - [60| . In the Gaussian model the relation holds 

r-TT 1 1/2 

(pT(a;))Gauss = ^(Pt (2;)) Gauss • (69) 

Of course, in no model considered so far such a factorized x- and transverse parton momentum dependence was ever 
observed, and in the bag model we do not observe it either. However, it is interesting to ask, for example, to which 
extent the relation is supported in a model. With {p'^{x)) defined in ([M]) we obtain for the expression on the RHS 
of the result plotted as dashed line in Fig. [5J:. The remarkable observation is that ((M)) is supported within an 
accuracy of 0(10 %) in the valence- a; region. We remark that this conclusion is insensitive to the way the integration 

constant in f'\^^'^ (x)reg is fixed, provided this is done in a reasonable way (see above). 

However, the bag model supports the Gaussian model much more than that, in the following sense. In the Gauss 
model we have ff{x,pT) = ff{x,0) exp(-p|,/(p|(x)) Gauss) where, by definition, /f(x,0) = /f (a;)/(7r(p|,(a;))Gauss)- 
When dealing with a model with non-Gauss-likc transverse momentum dependence, this can be used to 'fit' the 
Gaussian width 

(p|(a;))Gauss = ^45rr (^0) 

such that the Gaussian model is exact at px ~ 0. By continuity arguments the Gaussian model can be expected to 
be a good approximation to the exact model results also for > in some vicinity close to pT ~ 0. The question 




FIG. 3: (a) (py) of ff as function of x. Solid line: computed using the exact definition in Eq. (|66p . Dashed line: using the 
Gauss model relation, Eq. H70[l. (b) and (c) (a;,PT) vs. pr for selected values of x. The thin dotted lines are the respective 
Gauss model approximations obtained from the Gaussian widths from Fig. [3^. 
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is: how large is the pr-region in which the Gaussian model with the width ([70]) will provide a useful approximation 
to the pT-dependence of the unpolarized TMD? 

Let us first compare the results for the mean transverse momentum square {p\{x)) defined in Eq. (|66p and the 
expression from the Gauss model Gauss, Eq. (ffO]). It has to be noticed that these are a priori completely 

different quantities. It is therefore remarkable that the results agree so well, especially for valence x, see Fig. [3^. 

In order to see to which extent the bag model is compatible with a Gauss-like shape of the transverse momentum 
distributions, we plot the pT-dependcncc of fi {x,pT) for < pr < Mn for selected values of a; = 0.1, 0.2, 0.3 in 
Fig.[5j3, and x = 0.4, 0.5, 0.6 in Fig.[3J;. In Figs.[3)D and c, we also plot the respective Gaussian approximations (as thin- 
dotted lines). The result is astonishing: In the valence- a; region the exact curves and their Gaussian approximations 
agree excellently! 

Not visible in Figs. |3Jd and c is that the first worthwhile mentioning deviations from the Gaussian behaviour start 
to become apparent only at larger px > 0.5 GeV. The crucial difference is obviously in the large-pT asymptotics: 
fi{x,pT) ^ a{x,pT)/pT with |a(a;,pT)| < const. (The function a{x,pT) oscillates, with a period defined by the 
periods of the spherical Bessel functions in (|14p. around some value which depends on the TMD but not on x. For 
all TMDs the respective functions a{x,pT) are bound from above and below.) 

Of course, non-perturbative models aiming at an effective description of the nucleon properties at hadronic scale 
are not able to address the large-pT region, where one may apply perturbative QCD. Effective models can, however, 
provide valuable insights for transverse momenta up to the order of magnitude of the hadronic scale, i.e. for px < Mj^. 
In this py-region the bag model supports the concept of a Gaussian distribution of the transverse parton momenta in 
the case of f1{x,pT)- 

Let us now turn our attention to other TMDs, keeping the discussion shorter after the detailed investigation of f^. 
The TMDs g', hi, e"^, gx, which exist also as collinear parton distribution functions, have the same large-p^ 
behavior like /f , and consequently also have divergent (l)-momcnts (which can be regularized as in the case of ). In 
contrast to this gj^, hj^i,, Z^*^, /i^p' behave like a{x,pT)/pT large pt and have well-defined convergent (l)-moments. 
Finally hf^, 9l'^7 .9t'' behave like a{x,pT)/PT and at large p^have well-defined convergent (l)-momcnts, too. 

Is it possible to approximate the pT-dependence also of the other TMDs in a Gaussian model? The answer is yes. 
Fig. m shows the Gaussian widths (p|(a;))Gauss as defined in Eq. ^ for the aU TMDs. (Widths of hfj^, /^«, g^'' 
are not shown, because they correspond to those of respectively h^"^ , h^^ thanks to the relations (j4TH43|) . Also 
the widths of gj^"^ = — ft,^ arc not shown, because of the discarded relation gj^"^ = —kz/AIpf h^^ implies the same 
PT-behaviour of h^^, gj^"^, h3p.) 

We observe a modest x-depcndence of the various (p^(a;))Gauss, see Figs.|3^-c. Important for the widths of 5^, h\ 
is that they are not larger than that of /f in order to comply with positivity, which is of course the case, see Fig. 2^. 
However, the widths of other TMDs are not bound in this way by (p|^ (a;)) Gauss of fl- In fact, the (p^(a;))Gauss of 
the pretzelosity distribution h^J^ exceeds the width of , see Fig. HJ:. Taken literally this would imply a violation of 
positivity, but we have to keep in mind that the (pf^ (a;)) Gauss only approximate the true p^-behavior of TMDS in the 
model, and the exact model results always satisfy positivity, see Sec. lIIIFl and App. IB] 
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FIG. 4: The Gaussian widths as defined in Eq. (|70p vs. pr for various TMDs. 
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Fig. [5] shows the pr-dependence of the TMDs for selected values of a;, chosen to optimize the clarity of the plots. 
The thin dotted lines are the respective Gauss model approximations obtained from the Gaussian widths from Fig. |4l 
We observe in general a good agreement, including even gj^^'- which has a zero in x, see Fig. |4l. 

Some TMDs, most notably for gf and e'' in Fig.|4K andHJ;, have for fixed values of a; a zero in pr- One is tempted to 
suspect a model-artifact which can be traced back to the oscillatory behaviour of the Bessel functions in However, 
in a covariant parton model calculation it was observed that gi{x,pT) becomes negative for some values of x, px (83| . 



gi(x,PT) (a) h'\(x,pj) (b) giri^.Pr) (c) 




0.2 0.4 Px/GeV o 0.2 0.4 Pi/GeV o 0.2 0.4 Pj/GeV 




hl(x,pj) (g) hj''(x,pj) (h) gi"(x,pj) (i) 




0.2 0.4 Pt/GcV 0.2 0.4 Px/GeV 0.2 0.4 Px/GeV 



FIG. 5: The pT-dependence of various TMDs for selected values of x. The thin dotted lines are the respective Gauss model 
approximations obtained from the Gaussian widths from Fig. U) 
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This detail may deserve further attention. But the effect is in any case small, and in phenomenological apphcations 
the total result, for example, in SIDIS after convolution with a fragmentation function, is strongly dominated by 
contributions from the valence- a; region and pr S 0-4 GeV, where the TMDs are sizable. 

Thus, to draw an intermediate conclusion, in the valence-a: region the bag model strongly supports the Gaussian 
model with a weakly a;- independent Gaussian widths (pf^ (x)) Gauss- 

Thus, the bag model results also encourage to use the so convenient Gaussian model for TMDs also in future (until 
the data will teach us the opposite). But how to use in practice these model predictions? Indeed, one cannot use the 
bag model predictions literally, since, for example, the Gaussian width of from bag model underestimates what is 
needed in phenomenology, see Eqs. (|67l) an d (f68l) . 

However, the pj'-broadcning mechanism |ll| is in lowest order approximation polarization independent. Therefore, 
model predictions for ratios of widths can be expected to be more reliable and useful for practical applications. In 
view of the weak x-dcpcndcncc of the (py (a;))Gauss wc chose z^, = 0.3 as a typical valence-x value and summarize the 
results for the (p^ (a;)) Gauss, in units of the width of /f, in Table HI 

It is worth to compare to the results from the light-cone constituent quark model [8l|, [s^ . In that model the 
wave-functions fall off with pr sufficiently fast, such that (p^) defined in Eq. (|64p exists. For the in that model 
{Pt) — 0-080 GeV^ is close to our {pxixv))Ganss — 0.077 GeV^. For gf, h\ the results are comparably similar. However, 
in the case of 55*7?, h^l^, h^^ the light-cone constituent quark model yields smaller widths compared to the bag model, 
see Table n Twist-3 TMDs were not studied in [Hiil. 



TMD j 


(Pt) 
Ref. [86] 


(pr(Xi,)) Gauss 

(bag, here) 


f? 


1 


1 


9l 


0.74 


0.71 


hi 


0.79 


0.85 


Hit 1 "iL 


0.74 


0.95 


h^" 
"it 


0.63 


1.11 


e" 




0.68 


5t 




0.84 


hi 




1.01 






1.11 






0.94 



TABLE I: Average transverse momentum squares in T-even TMDs from light-cone constituent quark model [s^l, and the bag 
model (results obtained here). The (py) from ,8^ are defined according to (|64[). The bag model results for the Gaussian widths 
are defined according to (|70p and taken at the valence- 2; point a;„ = 0.3. All results are in units of the respective value for 
which is = 0.080 GeV^ in the case of and {pI{xv)}g^I^^ = 0.077 GeV^ in the case of the bag model. 
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C. WW- type approximations 

By exploring the equations of motion, twist-3 parton distributions typically can be decomposed into pieces related 
to leading-twist TMDs, current quark mass terms, and quark-gluon-quark correlators. The latter are often referred 
to as 'pure twist-3' or 'interaction dependent' terms, and are marked by a tilde. For T-even TMDs one obtains [l3l 
(we suppress the arguments x and fc^) 

(71) 
(72) 
(73) 
(74) 

IVl 

-•2 



xe'' 




M 








+ /f , 








+ .9? + 


M 


xgrp 


= xgj. 


1 ig 
+ .9iT - 


^ M 



t2 



X 



K = xhl^^^hil + ^gl, (76) 



-•2 



?2 



a: 



4' = xh^+hl^^htS, (78) 



2 



.,- = .^- + 5/.?-^^/^^^ (79) 

where in the last equation the notation is used g'^ = g^ — g^'^ and analog for g'r^ [2^ . If we systematically 
assume that pure twist-3 and quark mass terms are small, which we indicate symbolically and generically by ©(e), 
and integrate over transverse momenta, then we obtain the following Wandzura-Wilczek-type approximations 



xe'^{x) 


= 0{e), 


(80) 


xf^\x) 


= /f(x)+0(£), 


(81) 




= gUx)+0{e), 


(82) 


xgr'^ix) 


= gt^ix)+0{e), 


(83) 


xg^{x) 


= gt^'^'^{x)+0{e), 


(84) 


xhlix) 


= -2hfl'^'^{x)+Oie), 


(85) 


xh^(x) 


= -h\{x)-h^^^\x)+0{e), 


(86) 


xh^'^{x) 


- h\{x)^h^^^\x) + 0{e), 


(87) 


xg'rix) 


- 0{e). 


(88) 



What these approximations have in common with the classic Wandzura-Wilczek (WW) approximation is that pure 
twist-3 and current quark mass terms are neglected. However, the neglected operators are different in all cases. 

Figs. |6^-[6l compare respectively the left- (solid lines) and right- hand-sides (dashed lines) of the Eqs. ([80H88)) . 
assuming that 0{e) is zero in each case. In all Figures the same scale is used in order to better compare the 
magnitudes of the various functions. First we observe that xe'^{x) and xg'^{x) are not zero, as one would expect on 
the basis on the WW- type approximations, though the functions are not large, see Figs.[S^ and[5J. 

Next we remark that the WW- type approximations for xf'^'^{x), a;/i' (x), xh^'^{x) can be considered as roughly 
supported by the bag model, sec Figs. |B)d, [Bf, Hh- In the remaining cases, however, the WW-type approximations 
work only modestly, for example Fig. [BJ;, or not at all, see Figs. IBt, [6ji, [6^. 

At first glance one might be surprised that the WW-type approximations are not exactly fulfilled in a no-gluon 
model with massless quarks, as apparently there are no contributions from quark-gluon and mass terms. However, 
here one has to recall that the WW-type approximations originate from applying QCD equations of motion, and 
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separating leading and subleading twist terms. In principle, one could repeat this game in the bag model, too. As the 
quarks are not free (but confined by the bag), one would consequently encounter certain "interaction dependent^^ tilde- 
terms in the model, too. These bag- model tilde-terms can be seen directly in the case of e'^{x), g'^ix) in Figs. [6^, 13. 
In the other plots in Fig. [6] they are apparent as the differences between the solid and dashed lines. It has been 
argued that the bag, which is a model for confinement, in some sense mimics the effects of gluons [l03l |. However, to 
which extent the bag-model interaction-dependent terms arc able to estimate reliably the QCD interaction-dependent 
terms, remains to be seen. 

In any case, it is interesting to observe that the bag model roughly supports the WW-type approximation for 
xf-^"{x), see Fig. |6}d, which played an important role in the interpretation [53 of the EMC data on the azimuthal 
asymmetry A™^"^ in unpolarizcd SIDIS [s^ as being due to the Cahn effect Q. 
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FIG. 6: Test of WW-type approximations in the bag model at the low scale. The solid lines show x times twist-3 TMDs, as 
given on the left-hand-sides on the Eqs. (|8UH88p . The dashed lines are twist-2 TMDs (or their transverse moments or linear 
combinations thereof) as given on the right-hand-sides of the Eqs. (|80H88(I . If the tilde- and mass-terms defined in terms of the 
QCD equations of motion vanished, the solid and dashed lines would coincide, see text. In Figs, (a) and (i) no dashed lines 
appear because the right-hand-sides of the Eqs. (|80|) and (|88p are in the WW-type approximation. 
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Having discussed the WW-type approximations for TMDs, whose usefuhiess rem ains to be tested it is interest- 
ing to have a look back on the "classic WW-approximations" for g^{x) and h'j^{x) |l03l . ll09j . which are distinguished 
from the WW-type approximations in t hat i i i the ir derivations the notion an d cornplic ations of transverse parton 
momenta does not need to be involved jl03l Il09l | , though can be considered jlOll Il24j | . These are therefore in a 
certain sense "collinear" approximations. They are given by 

sU^)"^^ f^9Uy), (89) 

Jx V 

hli.)"^^ 2.r^hUy). (90) 
J X y 

Figs. [7^ and [7}3 show to which extent the WW approximations are supported by the bag model: moderately in the 
case of g!^{x), and somewhat better in the case of /i' (.t). (In Figs. [7^ andlTt) we compare x times the functions and 
their WW-approximations, because at small x < 0.1 the approximations (|89 1 190| are poorly supported, which should 
not worry us as the bag model is expected to be more meaningful in the valence- a; region, see Sec.JlTI) 

Finall y, we remark that by combining the classic WW approximations for g^{x) in Eq. (|89p [103] . and h'j^{x) in 
Eq. (|90| jl03l |. with respectively the WW-type approximations in Eqs. ([84)) and (f55|) . one obtains in principle a further 
class of approximations [9^, which we shall denote by "WW & type" (short for WW and WW-type) approximations. 
These approximations relate leading twist TMDs to twist-2 parton distributions as follows [H, 

xf'^gliy), (91) 

J X y 

x^f^Kiy). (92) 

J X y 

These approximations were used in literature in order to make estimates for certain double (97j and single (98j spin 
asymmetries in SIDIS. In Figs. [71: and[7Ji we test the quality of these approximations in the bag model at the low 
scale. In both cases we observe that the approximations tend to overestimate the magnitude of the true model results 
for g^,^^°'{x) and hfj^^'^°'{x), somewhat more in the former case and less in the latter case. 

It is, however, difficult to suspect on the basis of these observations at the low scale of the bag model, whether the 
predictions from Refs. [13, [sll will similarly overestimate data. The evolution effects may play an important role. 
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FIG. 7: (a) and (b): The test of the classic WW approximations for g^{x) QO^ and hl{x) [10a |. see Eqs. ((2111501), in bag model, 
(c) and (d): Test of the approximations for g^^^^''{x) and h^j^^'^{x) in Eqs. (|91l I92|) . which result from combining Eqs. (|84II85|) 
and (|89l I90p . AH results refer to the low bag model scale. 
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Though the precise evolution pattern of polarized transverse momentum dependent distribution functions is not yet 
fully understood, it seems reasonable to expect the evolution from a low hadronic to experimentally relevant scales will 
"shift" the a;-shape of the TMDs towards smaller x, see also [1^. In other words, the region of valence- and large- a; in 
the bag model at its low scale, say 0.3 < a; < 0.8, could be what will be just relevant at higher scales in the experiments 
at COMPASS, Jefferson Lab or HERMES. Remarkably, in the region of larger x > 0.3 the approximations ([HTl 
work reasonably well, sec Figs. [71; and[7li. 



D. Orbital angular momentum 



For completeness we include here also the bag model results for the nucleon spin decomposition. Of course, it is well 
known that in the bag model 65 % of t he nucleon spin is due to the intrinsic spin of the quarks, and the remaining 35% 
are due to quark orbital moti on I l04ll . see also the study of this issue in the bag model in the context of generalized 
parton distribution functions [l25|. However, what is new here is that we obtain this information from the prctzelosity 
distribution, namely 

2LI = -2 J dx hf^'^^^'ix) = 0.46 , 2L^ = -2 / dx hf^'^^'^ix) = -0.11 , 2L^ = 2L^ + 2LI = 0.35 , (93) 

2Sl = J dx g'^ix) = 0.87 , 2S'3 = / dx gf{x) = -0.22 , 2S'| ^ 2Sl + 2Sl = 0.65 , (94) 

2Jl - 2Ll + 2Sl = ^ , 2Jl =2Ll + 2Sl^-^, 23% - 2L\ + 2^| = 1 . (95) 

This is a typical result for relativistic quark models at low hadronic scales |l26l . Il27j | . We remark that the MIT bag 
also reasonably well describes the axial coupling constant: the Bjorken sum rule yields gA = J dx [g^ — gi){x) = 1.09 
vs. 1.26 in experiment. (At this point it is in teresting t o no te that experimentally the Bjorken sum rule is verified 
with not much better accuracy than that |128| . see also |129| for a recent compariso n of QCD and data.) 

However, in order to comp are s uch numbers to phenomenology or lattice QCD |l30| one needs to carefully take 



into account evolution effects |13l| which, in the case of bag model, is supposed to start at a very low hadronic scal e 



[132| . Evolution techniques possibly suitable for that were discussed in the context of the Bjorken sum rule in |12£ 
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VI. CONCLUSIONS 



We presented a study of leading- and subleading-twist TMDs in the MIT bag model. Since this model laeks explicit 
gluon degrees of freedom, the Wilson-link needed in QCD to render the definitions of TMDs color gauge invariant 
is absent. As a consequence T-odd TMDs vanish. Attempts to model T-odd TMDs in the bag were presented in 
(oH - losj . In this work we have focused on the 14 T-even TMDs (6 leading- and 8 subleading-twist). 

Another consequence of the absence of the Wilson- link (in an y effectiv e approach with global color symmetry only) , 
is that certain relations hold among TMDs, the so-called LIRs IHlio^. There are 5 such LIRs among the 14 T-even 
leading- and subleading-twist TMDs, and we have proven that they are all satisfied in the bag model. 

Recently further relations among TMDs were found in models. One of the motivations of this work was therefore 
to shed some light on how such quark model relations arise. We have shown that in the MIT bag model there are not 
more and not less than 9 linear and 2 non-linear relations among the 14 T-even leading- and subleading-twist TMDs. 
We reviewed in detail that some of these relations are supported also in other quark models [zl, Ull, [s^, . 

One of those linear relations, found in our previous bag model stu dy [soj . connects the difference of gf and hf to the 
(l)-moment of pretzelosity. It was confirmed in several other [Sllll^,!!^ though not all [s^ quark models. What makes 
this relation particularly interesting, is the observation (s^ that such a difference between helicity and trans versify 



distributions is related in a light-cone SU(6) quark-diquark model |133| to quark orbital angular momentum (0AM). 

Although intuitively the idea of quark orbital motion is associated with TMDs, this is to best of our knowledge the 
first 'rigorous' connection of a TMD to 0AM — in a model, of course. While in gauge theories there is no gauge- 
invariant definition of an 0AM field operator, in quark models the situation is simpler and there are no ambiguities. 
Another important motivation for our study was therefore to investigate whether there exists any connection between 
pretzelosity and quark 0AM in the bag model. The answer is yes, agrees with the findings of [8J1, and reads 

Ll = (-1) J dx hf^^^^ix) . (96) 

Thus, by measuring the single spin asymmetry A^^^"^ due to hf^ [l3 one could access information on 0AM. Of 
course, the relation ([M)) is model-dependent. But it is supported by two independent approaches, bag model (here) 



and light-cone SU(6) quark-diquark model [8J|. Moreover, at least in the context of bag model, the information on 
0AM gai ned f rom pretzelosity, Eq. (|96|) . is equivalent to what one can learn from generalized parton distribution 
functions [l25| . It will be exciting to see to what extent experimental information on TMDs and generalized parton 
distribution will, on the basis of a quark model interpretation, converge to give a compatible picture of 0AM. 

The pretzelosity distribution /i^^ seems to play in this context a central role. It is also related to the non-sphericity 
of the spin-distribution in the nucleon [s^ . It is interesting to ask, whether a quark model relation of the type (j96p 
may inspire a way to establish a rigorous connection between TMDs and 0AM in QCD. The task is demanding, as 
we observe (|96p on the level of matrix elements only. Further studies in effective approaches, and numerical results 
from lattice QCD could provide valuable insights. For first attempts to study TMDs on a lattice see [s^]- 

The third main result of this work concerns practical aspects which are of interest for phenomenology. For example, 
in many phenomenological studies it is assumed that in TMDs the x- and pT^-dependencies factorize, and the latter is 
"Gaussian" . Many authors have stressed that in their models the x- and pr-dependencies of TMDs are non-factorizing 
and non-Gaussian. But in practical applications the Gaussian Ansatz works with a useful accuracy, e.g. [58l - l60| . 

How to reconcile these observations? In order to address this question, we introduced the notion of a (in general 
x-dependent) 'Gaussian width' which can be applied to any TMD. This effective Gaussian width is designed such that 
it reproduces the px-dependence of the TMD exactly in a vicinity of pt = by definition. Although in the bag we 
also observe non-factorizing, non-Gaussian x- and pT-dependencies, in this way we made two interesting observations. 
In the valence- a; region (and we speak here about a low hadronic scale), this effective 'Gaussian width' turns out to 
be only weakly x-dependent. Moreover, such an effective Gaussian Ansatz approximates reasonably well the exact 
model results not only in the vicinity of px = 0, but up to pt ^ O(Mjv). 

This is good news for phenomenological studies for two reasons: azimuthal asymmetries in Drell-Yan or SIDIS are 
dominated by intrinsic transverse parton momenta (sl-fioj. and one expects azimuthal and (single) spin phenomena 
to be sizable in the valence- a; region. And this is where we find the Gaussian Ansatz to be a useful approximation. 
Surely, care is required when sea-quark effects start to play a role, and for a precision treatment of transver se m omenta 



one has to use a rigorous approach such as the CoUins-Soper-Sterman formalism [11[ as i mplem ented in |134{ . 

Finally, we used the model results to test the Wandzura-Wilczek-type approximations |97l - llOCl{ which were suggested 
as, at the presently early stage, useful tools for first interpretations of data. These approximations consist in neglect 
'pure- interaction-dependent' terms (in QCD: quark-gluon correlations, in the bag model: bag boundary conditions). 
We observe that, for some TMDs these are fair approximations in the valence-x region. 

To conclude, though obtained in a simple model, our results bear many interesting insights, and we hope they will 
stimulate further studies in quark models. 
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Appendix A: Proofs of LIRs 



In this Appendix we prove the LIRs (p4l - |38l) . Thanks to the model relations discussed in Sec. IIIIBl we do not 
need to prove every LIR. For example, if we prove the LIR (p4| then also the LIR psp holds due to the relations 
(|42l l45l |46| , which can be seen conveniently by combining (|45| |46|) to form (j53p . Similarly, if we prove the LIR ((36|) 
then it is clear, that also the LIR (|571) holds, due to the model relations (^51 H^ . Finally, (|38p is evidently true, c.f. 
the model relation (|47)) . Thus, it is sufficient to demonstrate, for example, that the LIRs (p4|) and ((36|) hold. 

In order to prove p6p we rewrite the expression for hf^^\x) in a convenient for our purposes way. Recalling that 
fc is a function of k± and = xM^ — lu/Rq according to p^ . we may write 



,-L(i)g 



(x) 



PqA 



2Ml 



^2Mltl{k) 



PqA 



dr 
2ti 



siT?et\{q) 



(Al) 



where we replaced the transverse momentum integral by an integration over the independent variables q with q = |(f |, 
using a delta-function / dr/(27r) e*^'^^^'?^)'^ = — fc~), and the spherical coordinates qz = qcosO and q\ = q^ sir? 9. 
Next we differentiate (jAl[) with respect to x, recalling that the x-dependence is 'hidden' only in kz according to (fT8|) 



with j^kz = Mff- We obtain 



dx 



hf^'^'ix) ^ P,A I d'q 



dr 
2^ 
dr 
2^ 



^i{kz, —q cos 9)t 



[iMNT) 



PqA J d'q 

PgA j d\ 

PqA j d^q 6{qz-kz) 



N 



^i{kz —q cos 9)i 



q A cos ( 



dr Mn_ ^i{k,-q cos e)T 
27r q 



-2cos0tl{q) 



2MNq, 



PqA J d^fc^ 

— /i^(x) 



2MNkz tl{k) 



sin^etliq) 



(A2) 



where we interchanged the order of the integrations and differentiation, which is legitimate in our case, integrated by 
parts in the third step, and finally recovered the expression for (— l)ft,^(x), which completes the proof of the LIR ([36]). 
The proof of the LIR is straightforward, and consists of repeating steps analog to (|A11 IA2[) . 
It is interesting to remark that a different version of the LIR (p6)) is the following 

/.«'(x) -i^^^-^^^'W, (A3) 

with the subtlety that the (2)-moment of prctzelosity is divergent, but its derivative with respect to x is finite, i.e. 
the same careful treatment is required as in the case of the (l)-moment of /f discussed in detail in Sec. IV Bl The 
numerical bag model results for TMDs satisfy all LIRs including the version (|A3[) . 
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Appendix B: Proofs of inequalities 



This Appendix contains the exphcit demonstrations that the quark TMDs from the bag model satisfy the inequahties 
([51H55|) . We rccaU that, in the present version of the model, the T-odd TMDs /-[^ and h^"^ are absent and that the 
inequalities for antiquarks are violated, see Sec.|lT]for a detailed discussion. 

In order to check the inequalities in ([M)) we work directly with the model expressions. With fc;, = fc^/^/s^ +fc2 we 
have -1 < A:^ < 1 which we use below in ((BT|) . -1 < (2fc^ - 1) < 1 we use in (|B2)) . and fc^ < 1 used in (|B3|) . 



ff{x,k^)=NgA 
gl{x,k^)^P,A 



+ 2kz toil + t\ 



> N,Aito-hf > 



tl + 2k^ toh + {2kl - l)tl 



tl + 2k^ toil + tf 



h1{x, k±) 



4 + 2k, toh + kl tl 



<PqA 



tl + 2k, toil + tl 



f^{x,kj_) , 



(Bl) 
(B2) 
(B3) 



The inequalities (|B21IB3p imply Igfix, k±)\ < ff{x, k±) and (a;, k±)\ < fl{x, k±), because for the nucleon in SU(6) 



Pr. 



Na 



< 1 



u, d. 



(B4) 



In order to check the inequalities ([55l [55)1 it is convenient to explore the model relations ([5^ - H7|) . For example, 
because of (|B4p the relation (p9|) immediately i mpli es that for the nucleon the SofFer inequality ((55)) is a true inequality 
\hi{x, k±)\ < k±)+gKx, fc^)), see also jll6l |. Similarly, by using the relations (f39ll45l) to eliminate transversity, 

we conclude with (jB4p that ([55)1 is a true inequality, i.e. \hf'^{x, fcj_)| < k±) — g\{x, fcj^)). 

In order to verify the inequalities (j571 [58| we explore the linear (p9l447|) and non-linear ([STl [52]) model relations. 
It is sufficient, thanks to the absence of the T-odd TMDs and h\^'^ and the relation (|35|), to prove one of these 



inequalities, let us say 

(1)9 



For that we multiply ([ST]) by (fc^/2M^)2, and eliminate transversity by means of the 



relation ([39]) and h-^j, [x, k±) by means of the relations ([39l 145]) . This yields 



(B5) 



where in the last step we explored (|B4[) . 

It is interesting to remark that, except for f1{x, k±) > 0, all the other inequalities in ([MH55)) are true inequalities 
for the proton in SU(6), i.e. they are never saturated. For oth er ba ryons in SU(6) saturations may occur, for example, 
for strange quarks in the A hyperon, where Ng = Ps — I, see [ll6j l for a related discussion. 
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